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Within the framework of macroscopic quantum electrodynamics in linear, causal media, we study 
the van der Waals potentials of ground-state atoms in planar magnetodielectric host media. Our 
investigation extends earlier ones in two aspects: It allows for the atom to be embedded in a medium, 
thus covers many more realistic systems; and it takes account of the local-field correction. Two- and 
three-layer configurations are treated in detail both analytically and numerically. It is shown that 
an interplay of electric and magnetic properties in neighbouring media may give rise to potential 
wells or walls. Local-field corrections as high as 80% are found. By calculating the full potential 
including the translationally invariant and variant parts, we propose a way to estimate the (finite) 
value of the dispersion potential at the surface between two media. Connection with earlier work 
intended for biological applications is established. 

PACS numbers: 34.35.+a, 12.20.-m, 42.50.Nn 42.50.Wk, 



I. INTRODUCTION 



The dispersive interaction between small neutral, un- 
polariscd particles (atoms, molecules etc.) and an un- 
charged macroscopic object is a manifestation of the 
quantum nature of the electromagnetic field and related 
to the zero-point energy |l| . In the first quantum electro- 
dynamical treatment of this van der Waals (vdW) type 
force , and in many other related works, the atom was 
assumed to be located in free space. In reality, the atom 
may be submerged in or be a part of a gas, liquid, or 
even solid — a situation typical in diverse fields such 
as colloid science 0, 0, surface engineering [gj], and 
biology 7]. As a particular example from biology, one 
can cite the transfer of a small particle diluted in cell 
plasma through a cell membrane caused by vdW forces 
[8J. A question that may arise when the atom is embed- 
ded in a medium is how the vdW interaction is to be cor- 
rected due to the difference between the local fields act- 
ing on the atom and the macroscopic fields averaged over 
a region which contains a great number of the medium 
constituents, thus ignoring the gaps between them. In 
Ref. Q, the vdW interaction between two ground-state 
atoms embedded in adjacent semi-infinite magnetodielec- 
tric media has been considered. Comparison of this result 
with that deduced from the Casimir force on a thin com- 
posite slab in front of a composite semi-infinite medium, 
both obeying the Clausius-Mossotti relation, suggests a 
hint on how to account for the local-field correction. This 
confirmed shortly after in Ref. [lCj on the basis on a 
macroscopic quantum electrodynamics theory and the 
real-cavity model [ll|. A general formula for the vdW 



potential in an arbitrary geometry has been derived in 
the form of a sum of a translationally invariant term and 
and a term containing the uncorrected scattering Green 
tensor modified by a local-field correction factor [1, [Toll . 
This general formula is applicable for meta-materials [I2j 
and will serve as a starting point for our treatment of the 
local-field corrected vdW interaction in a stratified mag- 
netodielectric. Note that a generalisation of the formal- 
ism has recently been used to study local-field corrected 
interactions of an excited atom with a ground-state one 
across an interface between two media [13J. 

The interaction between a neutral atom and a mate- 
rial surface is customarily split into two parts: a short- 
range repulsive part, significant when the atomic valence 
electrons overlap with the surface, and a longer-range 
dispersive vdW part 0, EE El El EI EI E|- The- 
ories that focus on the vdW interaction commonly ne- 
glect the first, thus are incapable of correctly predicting 
the behaviour of the interaction potential at very short 
distances. They typically yield a power law of z~^ (za, 
atom-surface distance) for materials with dominant elec- 
tric properties, and a power law of z^ 1 for purely mag- 
netic materials, which lead to divergent values for the 
vdW potential right on the surface d [2l|, [22j . Due to 
the importance of phenomena such as physisorption and 
transport of particles through interfaces and membranes, 
much effort has been devoted to a better treatment of the 
atom-surface interaction when the two are at very short 
distances. With respect to physisorption, this has been 
done via introducing a reference plane [23] ; via charac- 
terising the material surface by a more realistic response 
function which includes spatial dispersion (24|, smooth 
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variation of the dielectric properties at the interface, and 
the contribution from d electrons to the screening [25j ; 
and via using an atom ic p olarisability going beyond the 
dipole approximation [l7], [H|. These studies typically 
produce finite values for the interaction potential at the 
interface (23J. 

In the present paper, the local-field corrected vdW in- 
teraction of a ground-state atom embedded in a planar, 
dispersing and absorbing magnetoelectric host medium 
is studied within the framework of macroscopic quantum 
electrodynamics and the real-cavity model. We propose 
a procedure that would allow for obtaining an estimate of 
the finite vdW energy of an atom right on the interface 
even within our predominantly macroscopic framework 
treating the atom within electric dipole approximation. 
The paper is organised as follows. In Sec. |TT] the main 
results concerning the local-field corrected vdW interac- 
tion of a ground-state atom and an arbitrary absorbing 
and dispersing macroscopic body are reviewed, and then 
applied to planar multilayer systems. The two-layer case 
is considered and a procedure for an estimation of the 
vdW potential at the interface is given in Sec. IIIII while 
the three-layer case is discussed in Sec. IIVI Connection 
with earlier work which focus on biological systems and 
concluding remarks are given in Sec. |V] 



II. THE MODEL 

We set the stage by briefly recalling the results for 
the local-field corrected vdW potential of a ground-state 
atom within an arbitrary geometry and applying them 
to a planar multilayer system. 



A. Local-field corrected vdW potentials in 
arbitrary geometries 

We consider a ground-state guest atom A located at 
ta in an absorbing and dispersing host medium of arbi- 
trary size and shape characterised by macroscopic er(r, u) 
and /j,(r,oj). To account for the local-field correction we 
employ the real-cavity model by assuming the atom to 
be surrounded by a small spherical free-space cavity of 
radius R c . The radius of the cavity is a measure of the 
distance between the guest atom and the surrounding 
host atoms 10]. By construction, the macroscopic quan- 
tities s(r,uj) and ^(r,cj) do not vary appreciably on the 
microscopic length scale R c . To apply the real cavity 
model one should keep in mind that ^/e(0)^(0)i? c should 
be small compared to the maximum of all characteristic 
atomic and medium wavelengths as well as to the separa- 
tion from any surface of the host medium. In particular, 
the application of model to metals is very problematic. 
Using second-order perturbation theory, the vdW poten- 
tial for such an atom can be written in the form II Oil 



where U\(ya) is constant throughout any homogenous 
region and accounts for all scattering processes within 
the cavity, 



U^va) 



4n 2 c 



(2) 



where 



CUM 





zh ( ^\z) 


' -e A (<d)h?\z) 




e A (w)h^\z) [zoji(zo)]' -ji(zo) 




1 



(3) 



[z =uiR c /c, z = n(uj)z Q , n(u) = y/e A (u)ii A {u), e A {u) = 

e(r A ,u), (J-a(u) = h(t a ,u)], with j x {x) and h^\x), re- 
spectively, being the first spherical Bessel function and 
the first spherical Hankel function of the first kind, 
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and a(u>) denoting the isotropic polarisability of the guest 
atom in lowest non-vanishing order of perturbation the- 
ory H|, 
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with cow — (Ef. — Eq)/% being the (unperturbed) atomic 
transition frequencies and d' fe = (fc|d|Z) being the atomic 
electric-dipole transition matrix elements. 

The term U2{ya) involves all interactions associated 
with the particular shape and size of the magnetodielec- 
tric host medium, 

U 2 (v A ) = ^-J°° da 2 £>i(W0Tr G« (v A ,r A , *0 
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where 
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The Green tensor G*- 1 -* (r, r', uS) of the electromagnetic 
field accounts for scattering at the inhomogeneities of 
the (unperturbed) magnetoelectric host medium while 
the factor D A {u>) comes from the local-field correction. 

Within the real-cavity model, the potentials U\{r A ) 
and U2{ya) are well approximated by their asymptotic 
limit of small cavity radii, i.e., we keep only the lead- 
ing non- vanishing order in y^e(0)/i(0)[w roax -R c /c where 
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Umax represents the maximum of characteristic atomic 
and medium frequencies (for details, see Ref. fiolp. 
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[j > 0, for j = set do = 0] with the abbreviation 
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In Eq. (TIB"]) . p(s) denotes p(s) polarisations. The reflec- 
tion coefficients obey the recursion relations 
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and 
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Here and in the following the dependence of e, [i and a 
on the z£ is suppressed for brevity. Note that in leading 
order the local- field factor [3e A /(2s A + l)] 2 depends on 
dielectric properties only. The associated (conservative) 
vdW force is given by 
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The cavity-induced part of the potential, U\ according 
to Eq. @ does not lead to a force but to an energy shift 
whose influence on the overall potential will be studied 
in Sees. HI] and HVl 



B. Atom in a magnetoelectric multi- layer system 

So far we have not been specific about the geometry 
of the macroscopic body. Consider a stack of n layers 
labeled by I (I = 1, . . . , n) of thicknesses di with planar 
parallel boundary surfaces, where e(r, u>) — £i{u>) and 
fi(r,u>) = Hi{ui) in layer I. The coordinate system is 
chosen such that the layers are perpendicular to the z 
axis and extend from z = to z = di for I ^ l,n and 
from z — to z = — oo (oo) for I = 1 (n). The scattering 
part of the Green tensor at imaginary frequencies for r 
and r' in layer j is given by (see, e.g., Ref. (22J) 

G (1) (r,r',i£) = J dV <q -<'- r '>G (1) (q,*,* / ,if) (12) 
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= 0], where the modulus of wave vector in the 
2 direction is given by 

A = v^ 2 + ? 2 ( 17 ) 

with fc;, which always appears in the form of kf in the 
Green tensor, Eq. (fT3|) . being 
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The s- and p-polarisation unit vectors are defined as 

1 
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[e q = q/q,q = |q|]. Substitution of Eqs. (fT2"j) and |T 
into Eq. {TOJ) leads to 
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where we have used the relations 
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to calculate the trace. It is worth noting that the term in 
curly brackets in Eq. (|20|) describes processes that involve 
an odd number of reflections at the interfaces while the 
second term accounts for an even number of reflections, 
as can be seen from 
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Obviously, the expression presented in Eq. (|25[> . which 
corresponds to even numbers of reflections, is indepen- 
dent of the position of the atom and rapidly decreases 
with increasing distance between two plates. Hence, it 
cannot lead to a force but only to a layer-dependent en- 
ergy shift which is small compared to the other terms. It 
is completely absent in the two-layer case since only single 
reflections at the interface occur. It also vanishes for suf- 
ficiently dilute media, an expansion for small Xj — e j~^ 
and Q — fij — 1 showing no contribution to linear or- 
der, which reflects the fact that at least three medium- 
assisted reflection processes are involved. On the con- 
trary, Eqs. (|2"5|) and (j2"4"|) depend on the position of the 
atom, where waves propagating to the left and to the 
right of the atom have to be distinguished. 

It is worth noting that all formulas in this section are 
valid for arbitrary (passive) magnetoelectric media, in- 
cluding mctamaterials and in particular, lefthanded ma- 
terials with simultaneously negative real parts of e and /i, 
due to the fact that y/£j(ui)fj,j(ui) has no branching points 
in the upper half of the complex frequency plane. How- 
ever, left-handed material properties being only realized 
in finite frequency windows, they can not be expected 
to have a strong influence on ground-state dispersion po- 
tentials which depend on the medium response at all fre- 
quencies; alternatively, this can be seen from the fact that 
the potential is expressible in terms of permittivities and 
permeabilities taken at imaginary frequencies which are 
always positive [H, H3] ■ The situation can change when 
the atom is in an excited state and the atom-field interac- 
tion resonantly depends on the medium response at the 
atomic transition frequencies [H, Hi| . 



III. VDW POTENTIAL NEAR AN INTERFACE 

Let us first apply the theory to the simplest case of a 
single interface between two homogeneous semi-infinite 
magnetoelectric half spaces, where we first concentrate 
on the position-dependent part of the potential respon- 
sible for the vdW force and then interpret the (layer- 
dependent) constant part. 



A. Position-dependent part 

For a planar magnetoelectric two-layer system with the 
guest atom placed in, say, layer 2, substitution of = 
and in accordance with Eqs. (fTBj) in Eq. (j2"U)) leads to 
the following expression for the position-dependent part 
of the vdW potential 
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which differs from its atom-in-free-space counterpart by 

the local-field correction factor ( 2 <r2+i ) an< ^ £ 2 =/= 1 
and Ai 2 ^ 1 [H 



1. Retarded limit 

It is instructive to study the potential in the retarded 
limit 



c c 
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where and Wj f being the minima of all relevant 
atomic transition and medium resonance frequencies, 
respectively. Introducing a new integration variable 
v = c/3 2 /t; and replacing a(i£) ~ a(0), £1,2 — £ i,2(0), 
and fJ-i^{iS,) — A*i,2(0), the integration over £ can be per- 
formed and we obtain for the potential 
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It can be proven that <9C 4 /<9ei(0) < 0, dd/d^Q) > 0, 
and dC 4 /dfi 2 (0) < 0. 

To deduce some physics from Eq. (f2"9")h we consider 
some limiting cases. Assuming that the atom is located 
in free space ^ 2 (0) = £2(0) = 1, it can be shown that, for 
a purely electric halfspace 1, 

C 4 [M0) - 1,M0) = Ma(Q) = 1] < (30) 
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and, for a purely magnetic halfspace 1 with /ii(0) > 1, 

C 4 [ei(0) = l,na{0) = l,e 2 (0) = 1] > 0. (31) 

The inequality (|30p means the atom is attracted toward 
the electric halfspace — the case which is commonly 
treated in earlier literature. On the other hand, the pos- 
itivity of C4 in Eq. (|31|) means the atom is repelled from 
the magnetic half space. This, coupled with the signs of 
the derivatives given below Eq. (f2l*"[) . imply that electric 
properties tend to make the potential attractive while 
magnetic ones tend to make the potential repulsive. 

Now if the atom is embedded in a material halfspace, 
while the opposite halfspace is empty /ii(0) = £i(0) = 1, 
it can be shown that, for a purely electric material, 

C 4 [mi(0) = l,£i(0) = l,/i 2 (0) - 1] > 0, (32) 

the atom is repulsed from the interface, while for a purely 
magnetic material with /i 2 (0) > 1, 

C 4 [mi(0) = Mi(0) = l,e 2 (0) = 1] < 0, (33) 

the atom is attracted towards the interface. Since 
dCi/dfi2{0) < 0, if one starts from a purely electric ma- 
terial which is accompanied by a repulsive potential, then 
enhances the magnetic properties of the material by in- 
creasing ^2(0), one would obtain with an attractive po- 
tential in accordance with Eq. (f3"3"]l . 

Another particular case is when the magnetodielectric 
contrast between the contacting media is small 

Ei(0)=e 2 (0)+x(0), (34) 

Mi(0)=M2(0) + C(0) (35) 

[x(0) < £2(0), C(0) < Ma(0)], one can further treat the 
integrals analytically by keeping only terms linear in \ 
and £, 
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The u-integration is straightforward and we arrive at 
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This result generalises the one obtained in Ref. [22j to the 
case of an atom embedded in a medium, with local-field 
correction included. It is richer in content than its atom- 
in-free-space counterpart. For example, when x(0) = 
£(0) and the atom in free space £2(0) = /i 2 (0) = 1, the 
potential is attractive, whereas when the atom is in a 
medium of e 2 (0)//t2(0) > 23/7, a repulsive potential can 
be realized. 



2. Non-retarded limit 

The non-retarded limit corresponds to atom-surface 
distances za small compared with the typical wave- 
lengths of the medium and the atomic system, 
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[711,2(0) = y/si, 2(0)^1,2(0); 0/4 and uij^, maxima of the 
relevant atomic transition and medium resonance fre- 
quencies]. The conditions (|3"5|) and (|4"0"]) imply 
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where we have used that the ^-integration is practically 
limited to a region where £ < u>\ M . We substitute 



1 = \jP2 - £2/i2C 2 /c 2 , 
qdq = /? 2 d/? 2 , 



0i = \J / ^ 2 /c 2 (ei/ii - 82^2) + 0: 
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in Eq. (|26|) and, on recalling Eq. (|4lj) . perform a leading- 
order Taylor expansion in £ 2 /(c 2 /3f )(ei/ii ~ £2/^2)- After 
carrying out the /3 2 -integral, we arrive at 
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(48) 



and Eq. ([4*2")) has been used to set exp(— 2^/£ 2 M2^a/c) 
equal to one. By putting e 2 = fi 2 = 1, Eqs . (f47l) and ([48]) 
reduce to those for an atom in free space [22j |. 

As can be seen from Eqs. p6|) - (|47| . one can distinguish 
two regimes having different power laws. The first regime 
is where the two contacting media have unequal electric 
properties. Then the first term in Eq. (|46[) dominates 
and the power law is z~^ . The atom is pulled toward 
(repelled from) the interface if the medium it is located 
in has stronger (weaker) electric properties than that on 
the other side of the interface. 
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In the case of equal electric properties C3 = and 
U%(za) = Ci/z A with 
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d^ 2 a ( 
V2£2 + l 
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The atom experiences a force which points away from 
(toward) the interface if the magnetic properties of the 
medium the atom is situated in are weaker (stronger) 
than those of the medium on the other side of the inter- 
face. 

Note that the dependence of the directions of the forces 
on the difference in strength of the medium responses 
is opposite in the two cases of dominantly electric and 
purely magnetic media. In both cases the strength of 
the force increases with increasing difference between the 
electric and magnetic parameters of the contacting me- 
dia. These results are consistent with earlier ones ||, [3(| • 
A more detailed comparison of our results with those pre- 
sented in Ref. [8| will be made in the last section. Since 
the coefficient C3 of the leading-order term depends on 
electric properties only, the influence of electric prop- 
erties on the behaviour of the potential tends to dom- 
inate that of the magnetic properties, except for the case 
when the electric properties of the neighboring media are 
similar. At more moderate distances, competing effects 
of electric and magnetic properties may create potential 
walls or wells, as is also evident from the numerical re- 
sults below. 



3. Numerical results 

To study the local-field corrected potential at interme- 
diate distances and to elucidate the combined influence 
of electric and magnetic properties of the media, we cal- 
culate the position dependent part U2{za) in accordance 
with Eq. (|26[) numerically. A two-level atom of transi- 
tion frequency u\q is assumed and material electric and 
magnetic properties are described using single-resonance 
Drude-Lorentz-type permittivities and permeabilities 
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Figure [T] (a) illustrates the atom-surface-distance de- 
pendence of the Ui[za) potential. The interface is at 
za = 0, the medium 1 is on the left, while the medium 
2 is on the right. The parameters e%, fi%, /12 are fixed 
whereas £2 takes on three different values of the plasma 
frequency wp e 2- In all three cases E\ and £2 have same 
(transverse) resonance frequencies uiTe and damping con- 
stants 7 e , but different plasma frequencies (i.e., e± =/= £2) 
and the zT 1 term in Eq. (|46[) is negligibly small. In case 
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FIG. 1: (a) Position-dependent part of the vdW potential ex- 
perienced by a ground-state two-level atom in a magnetoelec- 
tric two-layer system as a function of atom-surface distance 
for fixed ei, fii, /12, and for wp^/^io = 1 (1), 0.4 (2), and 0.2 
(3). Solid lines denote the potentials with the local-field cor- 
rection, while dashed lines represent those without. Other pa- 
rameters are lotei/luio = UTc.2/^10 = 1.03, u>p e i/u)io = 0.75, 

WTml/WlO = UJTm2 /Wio = 1, U>p m l/^10 = 2.3, LOp m 2/uiW = 

0.4, 7 m i,2/wio = 7ei,2/^io = 0.001, and the cavity radius 
is R c uiio/c = 0.01. (b) Difference between local-field cor- 
rected and uncorrected (position-dependent) vdW potential 
AU2 versus atom-surface distance where the solid, dashed, 
and dotted lines refer to the curves (1), (2), and (3), respec- 
tively. 



(1) £2 > £i, it can be seen from the figure that the po- 
tential at very short atom-surface distances is repulsive 
in medium 2 and attractive in medium 1, in consistency 
with the analytical result (|4U|) and (|4"T)l . Similarly, in 
cases (2) and (3) £2 < £1, the potential is attractive in 
medium 2 while repulsive in medium 1. As the atom- 
surface distance increases, the second term — in the po- 
tential (|46| gradually comes into play and if the magnetic 
properties are strong enough, may switch the sign of the 
potential and create potential walls or wells in the pro- 
cess, as is clearly visible in case (2). 

To show how the net effect of the local-field correc- 
tion depends on the distance and on the properties of 
the media surrounding the guest atom, we plot the un- 
corrected potential by dashed lines in Fig.[T](a), and the 
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£2(0) 

FIG. 2: Position-dependent part of the vdW potential as 
a function of the static permittivity £2(0) (more specifi- 
cally wp e 2/wio) for two values of the atom-surface distance 
zau)io/c — 0.01 (scale to the left) and zau)io/c = 3 (scale to 
the right). Solid lines are with the local-field correction while 
dashed lines are without one. Other parameters are the same 
as in Fig.fj] Note the arrows which indicate the ordinate scale 
to be used. 



difference between the corrected and uncorrected results 
in Fig. [1] (b) . The ratio between the corrected and un- 
corrected results is not always a good measure of the dif- 
ference between the two because one of them can vanish. 
The local-field correction factor [3e,-/(2£j + l)] 2 is posi- 
tive, larger than one, and increases with Ej [j indicative 
of the layer containing the guest atom] . ft approaches the 
maximum value of 9/4 as Ej —+00. Note that a larger- 
than-unity local-field correction factor does not necessar- 
ily lead to an enhancement of the potential because the 
uncorrected factor in the integrand can change sign as 
£ varies. The local- field correction has a clear-cut ef- 
fect of increasing or decreasing the potential only when 
the uncorrected integrand is purely repulsive or attrac- 
tive, for which cases (1) and (3) can serve as examples. 
In the middle case (2), the two curves with and with- 
out local-field correction cross, that is, there exists an 
atom-surface distance at which the effect of the local- 
field correction is canceled out due to the ^-integration. 
In addition, one can notice that the local-field correc- 
tion leads a small shift of the position of the peak. Fig- 
ure Q] (b) which shows the difference between the local- 
field corrected and uncorrected potentials, reveals quite 
significant corrections of up to 30% of the uncorrected 
values. 

The behaviour of the local-field corrected vdW poten- 
tial with respect to the static permittivity of the medium 
the atom is embedded in is shown in Fig. [5] for two differ- 
ent values of the atom-surface distance. Within the scale 
of the figure, the curves for the larger distance from the 
interface peak at certain values of £2(0). The positions of 
the peaks are different due to the effects of the local field. 
As wp e 2/wio and as a consequence £2(0) increases, an in- 
spection of the figure reveals that the ratio between the 
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FIG. 3: The exact layer-dependent constant part of the po- 
tential (solid line), Eq. J2J, and approximate results (dashed 
line), Eq. are shown as functions of the real-cavity radius. 
The upper pair of curves shows ~Ui(za)12tv 2 soc 3 / (Lof Q \dio\ 2 ) 
(the sign has been reversed so that a logarithmic scale can be 
used) for a pure electric material with uip e 2/ui\Q = 0.4, while 
the lower pair of curves shows !7i(zA)127r 2 £oc 3 /( w io|dio| 2 ) for 
a pure magnetic material with uJpmi/wio = 0.4. All other pa- 
rameters are the same as in Fig. [T] The radius of the cavity 
R c lj-lo/c starts from 0.001. 



corrected and uncorrected curves tends to the static value 
of the local-field correction factor [3£ 2 (0)/(2e 2 (0) + l)] 2 
(which lies between 1 and 9/4), in agreement with the 
analytical analysis given in Sec. IIII A II For the smaller- 
value of the atom-surface distance zaloiq/c — 0.01, a 
crossing point between the corrected and uncorrected 
curves is observed, where the local-field correction pro- 
duces no net change. 



B. Layer-dependent constant part 

Fig.[3]shows the dependence of the constant part U\ of 
the potential on the real-cavity radius R c . To gain some 
insight, we consider the two limiting cases of a purely 
electric and a purely magnetic material. According to 
the analytic result @, in the first case the leading term 
is proportional to [i? c a->io/c] -3 , while in the second case 
the leading term is proportional to [R c u>iq/ 'c] -1 . Thus for 
small enough R c <j}\q/c, \U\ \ for a pure electric material is 
generally larger than that for a pure magnetic material 
- a fact that is confirmed by the figure. It can be seen 
that U\ < in the first case and U\ > in the second 
case, again in agreement with Eq. ([9]). The figure also 
indicates that the magnitude of U± , which is entirely due 
to the local-field correction, decreases with an increasing 
real-cavity radius, that is, the effects of the local field 
becomes weaker as the medium becomes more dilute. 

For comparison, we also plot the potential V\ accord- 
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ing to the approximate result ((H) (dashed curves). It 
can be seen that for the parameters used in the figure, 
the approximate result reproduces quite well the exact 
one, especially in the case of a pure electric material. 
The agreement in the case of a pure magnetic material 
is good for very small RqUio/c, but worsens as R c ujio/c 
increases. 



Total vdW potential and its value at the 
interface 
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We have separately investigated the position-depen- 
dent part U2{za) of the potential, which determines 
the force acting on the atom, and the constant, layer- 
dependent part Ui , which is related to the local- field cor- 
rection. For problems such as the transfer of an atom or 
a small particle through an interface, it is of relevance 
to evaluate the potential right at the interface. For this 
purpose, the total value of the potential is needed. In 
Fig. [3J we have calculated U\ + ^(^a) on both sides of 
the interface with the medium 1 fixed while the medium 
2 varying from vacuum to a more dense medium with bal- 
anced electric and magnetic properties. The case repre- 
sented by dashed line is nothing else rather than case (2) 
in Fig. [TJ Only very short distances \J\emi\zAUw /c«l 
are presented. Since £\ ^ £2 in general, the position 
dependence of the potential at short distances is mostly 
determined by the term which contains £1 — £2 in 

Z A 

the integrand. Numerical results in the figure are consis- 
tent with this and show that the potential is attractive 
(repulsive) if the atom is located in a medium which is 
electrically more dilute (dense) than that in the opposite 
side of the interface. Additional structures in U2 like po- 
tential wells or walls are typically overwhelmed by the 
magnitude of U\. Baring the visual suppression of the 
potential wells or walls in U2 due to large magnitude of 
Ui, a full potential is more straightforward than U2 alone 
in predicting the movement of an atom across a surface. 
Take, for example, case (2) in Fig.[TJ An atom located in 
layer 2 close to the surface will be attracted to it, and if 
the atom can cross the interface, it will be pushed further 
away from the surface into layer 1 . Fig. 2] (dashed line) 
provides us with some additional information. It shows 
explicitly that the total potential in layer 2 is higher that 
that in layer 1, by giving the difference between the po- 
tentials in the two layers; and it may also help to estimate 
the amount of energy required for the atom to penetrate 
into layer 1. 

In Fig. [U we have not displayed the results for dis- 
tances \za\ < R c \/\£fJ-\ where the real-cavity model can 
no longer be applied. This gives rise to a gap between 
the potentials on the two sides of the interface. To ex- 
tend our theory so that it can be employed to estimate 



FIG. 4: Local-field corrected total vdW potential of a ground- 
state two-level atom in a magnetodielectric two-layer sys- 
tem as a function of the atom-surface distance. Different 
curves are for different (equal) electric and magnetic coupling 
strengths of the medium 2: ujp m 2/uJio = ujp^/^io = (solid 
line), 0.4 (dashed line) and 1 (dotted line). Other parameters 
are the same as in Fig. [1] 



the potential right at the interface, we suggest that 



U(z A = 0) = ~ [U{R C ) + U(-R c )] 
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(52) 



[recall Eqs. ©, (|46]), and (J47J)]. Visually, this means first 
plotting the potential as a function of the atomic position 
up to distances \za\ — R c (Rc being the radius of the 
cavity in the real-cavity model) , then connecting the two 
loose ends on the two sides of the interface to find the 
value of the potential at za = 0. Our result is remarkably 
similar to what has been obtained by calculating the on- 
surface potential of a molecule of finite size s [3l| , 
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(53) 



The second terms in Eqs. (1521 and ([53^1 . which represent 
the interface contribution to the potential, agree when 
setting s = ({/l6/3/n- 1 / 6 )R c w 1AR C and neglecting 
the local-field correction in Eq. (|52"|) which was not con- 
sidered in Ref. [311 ]. The first terms, which represent 
bulk contributions from the two interfacing media, dif- 
fer in the two approaches, where Eq. (|53|) still contains 
self-energy contributions which do not vanish in the vac- 
uum case Si — I while Eq. (|52p does not. Our result 
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thus represents an improvement of the previous one (|53[) 
in that local-field corrections are taken into account and 
self-energy contributions have consistently been removed. 



IV. THREE-LAYER SYSTEM 

A system consisting of an atom in a three-layer planar 
structure can serve as a prototype for theproblem of a 
small particle near or inside a membrane [8J. The trans- 
lationally invariant part and position-dependent part of 
the vdW potential can again be determined in accordance 
with Eqs. ([2]) and ([20]) , respectively. If the atom is in one 
of the two outer layers, Eq. (|2"D)) simplifies greatly via 
either rj + = or rj_ =0. If the atom is in the middle 
layer, the overall form of the potential Ui{za) remains 
as in Eq. (120]) and there exists a term which contains 
the product rj_rj + and is position-independent. This 
position-independent term in Ui{za) is irrelevant when 
it is the force which is concerned, but has to be kept, 
together with Ui, if the potential at the surfaces is of in- 
terest. Since the formulas are complicated, we resort to 
numerical computation. 

Figure [5] shows the behaviour of the vdW potential for 
an atom placed in an asymmetric (left column) and a 
symmetric (right column) three-layer magnetodielectric 
structure for different thicknesses of the middle layer. 
Note that the parameters for layer 1 are the same as 
those for layer 1 in Fig. Q] and the parameters for layer 

2 are the same as those for layer 2, case (2), in Fig. [2 
That is, the interface 1-2 here is the same as the interface 
1-2 in Fig. [2(2). In the asymmetric configuration, layer 

3 is vacuum while in the symmetric configuration, layer 
3 and 1 have the same characteristics. 

Let us first analyse the asymmetric configuration and 
see how the presence of a third layer 3 on the right, which 
is vacuum, affects the behaviour of the potential near the 

1- 2 interface. In case (a) where the middle-layer thick- 
ness dzWio/c — 5 is the largest among three cases, the 
behaviour of the potential U2{za) in the boundary re- 
gions is similar to that in the two-layer systems. Near 
the 1-2 interface, we find a potential well in layer 1 and 
a potential wall in layer 2, just as in Fig. [2(2). Near the 

2- 3 interface the potential is repulsive in the more dense 
medium 2 and attractive in the more dilute medium 3. A 
new feature appears around the center of the middle layer 
where a finite potential wall on the left and an attractive 
one towards the right interface combine to a potential 
well. Clearly, if an atom at rest is put in the well, it will 
remain there. When the thickness of the middle layer 2 
is reduced, the well becomes more shallow and eventu- 
ally disappears, as is visible in Figs. E](b) an d (c). For 
the parameters used in the figure, the potential well in 
the middle layer occurs when d,20Jw/c 3> 1 and is over- 
whelmed when d^io/c ~ 1. The curves for the potential 
can help one to predict the movement of an atom located 
near or inside a membrane (layer 2). For example, let an 
atom be initially located in layer 3 (vacuum) of case (a). 
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FIG. 5: Position-dependent part of the vdW potential as a 
function of zaujio /c for three different thicknesses of the mid- 
dle layer d,2LU\o/c — 5 [(a) and (d)], 2 [(b) and (e)], and 1 [(c) 
and (f)]. For the left column, the configuration is asymmet- 
ric With LOTel/^10 = 1-03, LJp e l/tJl0 = 0.75, = 
1, LOpml/cOlo — 2.3, LOTm2/(^10 = 1, WPmZ/WlO = 0.4, 
LOTe2/iOlO = 1-03, U)p e 2l^W = 0.4, 7m2,l = 7e2,l = O.OOltJlO, 

£3 = M3 — Ij an d RcUw/c = 0.01. For the right column, the 
configuration is symmetric with the vacuum in layer 3 being 
replaced by a medium of the same characteristics as those 
of layer 1. The curves without the local- field correction are 
shown by dashed lines. 



First, it will be attracted to the 2-3 interface. If it can 
be transported through the interface and gather enough 
momentum going down the slope, it can pass the finite 
potential wall and is then attracted to the 1-2 interface. 
After this interface, and maybe some oscillations, it will 
be suspended in the potential well near the surface. 

We now turn to the symmetric configuration (right col- 
umn of Fig. [5]) where the middle layer 2 is sandwiched 
between two identical layers 1 and 3. The behaviour of 
the potential on the right is just a mirror image of that 
on the left with the mirror plane being the one paral- 
lel to the surface and passing through the centre of the 
middle layer. When the thickness of layer 2 is largest 
d,2^w/c = 5 [case (d)], we see in layer 2 a combination of 
two potential walls as found in Fig. [2 with a well in the 
middle. Thus even as the middle layer is in general less 
dense than the two surrounding layers in both electric 
and magnetic aspects, there exists a possibility that an 
atom initially located in the potential well remains there. 
With decreasing thickness c?2 > the bottom of the well rises 
and the two walls eventually merge into one [cases (e) and 
(f)]. Clearly, in the cases (e) and (f), any atom initially 
situated in the middle layer will be transported to the 
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neighboring layers. Numerical computation also shows 
that the magnitude of the position-independent term in 
U2(za), i-e., the last term in Eq. ([20|) . is negligible com- 
pared to the position-independent terms, due to the small 
exponential factor. 

In Fig. [5j the uncorrected potentials are plotted by 
dashed lines. In the two outer layers, the effects of the 
local field almost remain the same as in the two-layer 
configuration, which can be explained by that the pres- 
ence of a third layer is screened by the middle layer. For 
the middle layer, we have found that typically, the local- 
field correction is most significant around the centre of 
the layer, as is most manifest in case (f) where a cor- 
rection [(^corrected — [^uncorrected) /^uncorrected] of more 

than 80% is observed. A variation of the middle layer 
thickness does not affect much the strength of the local 
field correction near the surfaces. This can be understood 
as resulting from the fact that when an atom is located 
very close to a surface, it will tend to see only the nearest 
neighboring layer. 

If one wish to know the potentials right on the surfaces, 
one would have to evaluate the full potential U\ + U^a) 
for atom-surface distances large enough such that the 
macroscopic theory applies, then use it as an input to 
the procedure proposed in the previous section. Since 
the position-independent part U± does not depend on the 
layer thicknesses, at each surface in a more-than-two- 
layer system, the results will be closely analogous to those 
of the two-layer case. 



V. DISCUSSION AND SUMMARY 

Our results might be of interest in biological applica- 
tions such as the transfer of a small molecule through a 
membrane from one cell to another. Earlier theories have 
been developed to describe the vdW interaction between 
molecules or small particles and a solvent medium |8j . In 
particular, such particles were modeled by a small dielec- 
tric sphere of radius R s and (macroscopic) permittivity 
e s . It was found that the nonretarded dispersion poten- 
tial of such a sphere near the interface of two dielectric 
media (with the sphere being situated in a medium of 
permittivity £2, and £1 denoting the permittivity of the 
medium on the far side of the interface) is given by 



U(z s ) 



16it 2 Eqz3 J e 2 £1 + £2 



(54) 



Here, z s is the distance from the sphere to the surface 
and 



, . . . £ s (u>) — £2(1^) 

e s (u>) + 2e 2 (w) 



(55) 



is the excess or effective polarisability [32J, |33| of the dis- 
solved particle in the medium. This potential already 
accounts for the fact that the (macroscopic) particle has 
a finite volume and can hence only move by displacing 



an equal volume of solvent from its path (with associ- 
ated pressure forces being present); the excess polaris- 
ability and hence also the potential must vanish when a 
dissolved particle has the same properties as the solvent. 

Equation (|54[) is the macroscopic counterpart of our 
microscopic equations (|46|) — (f48|) . where the position- 
dependent part of our potential for purely dielectric sol- 
vent media reads 



U 2 (z A ) = - 



n 



16n 2 e z A j Q 



£2 \££2 



3£9 



£l - £2 



£\ 



£2 
(56) 

The results from a microscopic description of the parti- 
cle as a system of bound point charges, leading to our 
Eq. (|56p . are thus formally very similar to those from a 
macroscopic model of the particle as a dielectric sphere, 
Eq. (|54)l . The main difference is the fact that the micro- 
scopic polarisability ^ together with a local-field cor- 
rection factor appears in Eq. (|56p . while the macroscopic 
excess polarisability ([55]) which effectively accounts for 
pressure forces enters Eq. ([54| . Depending on the size of 
the immersed particle, one of the two models may provide 
a more realistic description: For very small particles like 
atoms or small molecules whose size is comparable to the 
free interspaces between the atoms forming the solvent 
medium, local-field effects are important while (macro- 
scopic) pressure forces can not even be defined, so that 
the microscopic result (|56"|) should be used. For larger 
molecules whose volume covers a region that would oth- 
erwise be occupied by a large number of solvent atoms, 
local-field effects can be neglected while pressure forces 
become relevant, so Eq. (|54[) should be given preference. 
As an additional difference, note that our microscopic 
calculation has also given layer-dependent constant con- 
tributions to the potential, which are absent in Eq. (|54[) . 

Similar considerations apply in the retarded limit, 
where the macroscopic potential of the small sphere was 
found to be [8[ 



U(z s ) 



23hc 



a s (0) ei(0)-e 2 (0) 



320tt 



W s 4 £ 3 2 /2 (0) ei(0)+e 2 (0) 



(57) 



whereas our microscopic result as given by Eqs. (j28f and 
(HHJ) reads 



U 2 (z A ) 



3hc 



«(0) 



64n 2 e z\ £ 



dy 
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3e 2 (0) 
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(58) 



[a = ei(0)/e2(0)]. In addition to the observations made 
for the nonretarded limit, the interface-dependent pro- 
portionality constants [i.e., the last factors in Eqs. |57|) 
and (|58|] are now also different in general, note that they 
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do agree in the limit of small dielectric contrast between 
the media, x(0) = £i(0) - e 2 (0) < e 2 (0). 

To conclude, we have the local-field corrected vdW po- 
tential of a ground state atom embedded in a planar mag- 
nctodielectric multilayer system by analytical and numer- 
ical means. The theory allows us to extend earlier studies 
of the same system in two important aspects: first, one 
can allow for the atom to be embedded in a medium thus 
making the theory applicable to a larger range of realistic 
situations; second, the effects of the local-field correction 
is elucidated. 

The formulas for the interaction potential have been 
derived for an arbitrary number of layers where the case 
of two- and three-layer systems have been studied in de- 
tail. We have shown that the potential can be decom- 
posed into two parts: A layer-dependent constant part 
which depends on the real-cavity radius, i.e., the density 
of the medium the atom is placed in; and a position- 
dependent part which contains the local-field correction 
as a factor in the integral over frequency. For the lat- 
ter, we have presented retarded and non-retarded limits 
and the considered case of neighboring media of similar 
properties. Distance laws have been reestablished with 
effects of the local field included where the local-field cor- 
rection has been found to be as high as 80% in certain 
cases. Further, numerical calculations show that an in- 
terplay between electric and magnetic properties of the 
neighboring media may lead to the appearance of poten- 
tial wells or walls near the surface. These structures are 
potentially helpful as a trapping mechanism. Although 
these structures are much less intense in magnitude than 
those occurring for an excited atom [28|, [29[, they are 
more permanent because an atom in the ground state 
has an infinitely long lifetime. 

The constant part of the potential, which originates 
from local-field effects, does not contribute to the vdW 



forces, but it can facilitate our understanding of the 
movement of an atom near an interface, and is instru- 
mental in our proposed estimate of the on-surface value 
of the potential: After calculating the total potential in 
both sides of an interface up to distances equal to the 
radius of the (real) cavity, beyond which a macroscopic 
model no longer applies, the average of the two poten- 
tial values taken at these distances can be regarded as 
the potential at the interface. Our procedure improves 
previous similar estimates by including local-field effects 
and consistently removing self-interactions. 

In the case of the three-layer system, emphasis was 
given to the influence of the thickness of the middle layer. 
While new features may arise for a very thin middle layer, 
the behaviour of the potential is simply those at the in- 
terfaces combined if the middle layer is thick enough. 

Our results, which have been obtained on the basis of 
an exact, microscopic model of the atom-field coupling, 
are complementary to previous more macroscopic dis- 
persion potentials of particles modelled as small spheres. 
The model of preference depends on the size of the parti- 
cle in the specific situations considered. In the future, ef- 
forts should be taken to find a dispersion potential which 
holds for all possible ranges of particle sizes and includes 
both the microscopic and macroscopic potentials as lim- 
iting cases. Our considerations may easily be extended 
to other geometries, such as spherically or cylindrically 
layered host media. 
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